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I Abstract. We prove the exponential decay of eigenfunctions of reductions of Brown- 
On' Ravenhall operators to arbitrary irreducible representations of rotation-reflection 
(~| , and permutation symmetry groups under the assumption that the corresponding 

' eigenvalues are below the essential spectrum. 



1. Introduction 



(N 
> 

^ ' The Brown-Ravenhall operator can be considered as the (multiparticle) Dirac operator 
lO ' projected to the positive spectral subspace of free particles. This operator was 
■ introduced in [1] as a Hamiltonian of quantum electrodynamics (QED) correct to the 
^ ' second order in the fine structure constant (see also [2j). The higher order corrections 

O ■ predicted by QED should thus be treated as perturbations. The Brown-Ravenhall model 
turns out to be a good candidate for this approach, as the recent rigorous results show. 
^ ' Indeed, it is bounded below even in the many-particle case for physically relevant nuclear 
c3 ■ charges [5], HI E], and the structure of its spectrum resembles the one of Schrodinger 
operator - the essential spectrum forms a semiaxis [SJ [TJ [H [9], possibly with some 
eigenvalues below ionization thresholds P, |9] . This is in a remarkable contrast to the 
many-particle Coulomb-Dirac operator which has essential spectrum on the whole real 
axis and no eigenvalues, but is sometimes used as a formal unperturbed Hamiltonian in 
some QED calculations. 

Having in mind the intention to consider the Brown-Ravenhall operator as an 
unperturbed intermediate model, it is very useful to have information on the rate 
of spatial decay of its eigenfunctions. In this article we prove that for systems of 
particles with electric charges of the same sign (we consider the potential energy 
of interactions with nuclei as external field) the eigenfunctions decay exponentially 
provided the corresponding eigenvalues are below the essential spectrum. This will also 
be proved for restrictions of the operator to subspaces of wavefunctions with certain 
rotation-reflection symmetries. 
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There are numerous results concerning the exponential decay of eigenfunctions of 
multiparticle Schrodinger operators, including anisotropic estimates and lower bounds. 
A very detailed analysis of the non-isotropic exponential decay of eigenfunctions of 
Schrodinger operators in terms of a metric in configuration space is presented in [10] . It 
is proved in [TT] that the upper bound of [10] is exact at least for the ground state. A 
very simple proof of the exponential decay, based on the approach of [10] can be found 
in [T2|, Lemma 6.2. 

As for relativistic operators, exponential decay of eigenfunctions is proved for 
one-particle Chandrasekhar operators [TSj [Hj and some projected multiparticle Dirac 
operators [15]. For one-particle Brown-Ravenhall atomic Hamiltonians the exponential 
decay of eigenfunctions was first obtained in [16] for coupling constants of Coulomb 
potential not exceeding In the recent preprint [17] the exponential decay of bigger 
rate is shown to hold pointwise for all one-electron atoms with subcritical or critical 
coupling constants. 

The paper is organized as follows. In Section [2] we introduce the Brown-Ravenhall 
model together with some auxiliary constructions and formulate the main result in 
Theorem 12.31 Then in Section [3] we discuss the relevant properties of the interaction 
potentials. The proof of Theorem 12.31 is presented in Section IH with the proofs of 
technical lemmata postponed until Sections [6]— [H In Section [5] we prepare these proofs 
recalling two useful theorems which give sufficient conditions for the boundedness of 
integral operators. Appendix contains a couple of properties of modified Bessel functions 
for reference. 

2. The model and the main result 

In the Hilbert space L2(M'^, C^) the Dirac operator describing a particle of mass m > 
is given by 

Dm = —ioL ■ V + /3m, 

where a. := (01,02,03) and (3 are the 4x4 Dirac matrices [T8|. The form domain of 
Dm is the Sobolev space if-'^/^(R^, C'^) and its spectrum is (—00, — m] U [m, +00). Let 
Km be the orthogonal projector onto the positive spectral subspace of Dm'- 

1 —ioL ■ V + I3m 
A^ := - H . 

2 2V-A + m2 

We consider a finite system of particles with positive masses m„, n = 1, . . . , A^. To 
simplify the notation we write D^ and A„ for Dm^ and A^^, and also for their tensor 
products with the identity operators in L2(M^,C^), e. g. 

n-l N n-1 N 

® / (g) Dmr, (g) ® / and ® / ® A„„ ® O /, 

j=l k=n+l j=l k=n+l 

respectively. 
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N 

Let f}N '■= ® A„L2(M'^, C^) be the Hilbert space with the inner product induced 

n=l 

N TV 

by the one on ® L2(R^C'^) = L2(R^^,C^ ). In this space the A^-particle Brown- 

n=l 

Ravenhall operator is formally defined by 



, N N ^ 

^T) — 1 Tl.-'l' 



(2.1) 



N 



N 

A 

n- 



with 

:= TT A„ = ® A 

n=l 
n=l 

Here the indices n and j indicate the particle on whose coordinates the corresponding 
operator acts. 

In (12. ip Vn and Unj are the operators of multiplication by the potential energy of 
interactions of the particles of the system with an external field and between themselves, 
respectively. In most applications to atomic and molecular physics Brown-Ravenhall 
operators are considered in the Born-Oppenheimer approximation. Then Vn is the 
potential energy of the n*^ particle in the electrostatic field of static nuclei 

K 



K(xn) :=e„ J] I (2.2) 



k=l 



where is the electric charge of the particle, and Zk and are the charges and positions 
of the nuclei. The interaction between the particles is given by the Coulomb potential 
energy 



t/„,(x„,x,):=^^^. (2.3) 



I -^J I 

We will assume that all the particles of the system have the same sign of electric charges 
e-n, n = 1, . . . , N, but otherwise they might be different, as happens for exotic atoms, 
where some electrons are replaced with muons or even hadrons. The spin of each particle 
is assumed to be equal to 1/2, as always with Dirac and Brown-Ravenhall operators. 
This implies that the particles of the system are fermions. According to the Pauli 
principle, if some of the particles are identical, the wavefunction of the system should 
be antisymmetric under their permutations. This means that the operator (12.11) should 
be restricted to the subspace of consisting of functions which transform according 
to a certain irreducible representation of a subgroup 11 of the symmetric group Sn 
generated by transpositions of identical particles. Let be the orthogonal projector in 
Sjn onto the space of such functions. We will denote the restriction of Hn to S^^ := P^S^ 
by H§. 

We will assume that the subcriticality condition 

min CnZk > -2(2/7r + 7r/2)"^ (2.4) 

n,fc 
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holds. According to [4J, 71^ (and thus 'H%) is bounded below even if we replace the 
strict inequality in (12. 4p by a non-strict. Violation of such non-strict inequality usually 
leads to the lack of boundedness below, as shown in [3j for the case of single nucleus. 
As far as Hn (or any of its restrictions) is bounded below, it can be defined via the 
corresponding quadratic form. 

It is convenient to reduce Tiff using the rotation-reflection symmetries of the 
system. Let 7 be an orthogonal transform in M^: the rotation around the axis directed 
along a unit vector through an angle (f^, possibly combined with the reflection 
X I— s> —X. The corresponding unitary operator acts on the functions ip G Sjn as (see 
[IE], Chapter 2) 

N 

(0^^)(xi, . . . , x^) = J] e-^'^-"--«"^(7-ixi, . . . , 7-ix^). 

n=l 

Here S„ = —\an A an is the spin operator acting on the spinor coordinates of the n^^ 
particle. The compact group of orthogonal transformations 7 such that commutes 
with Vn and Unj for all ra, j = 1, . . . , (and thus with Hff) we denote by T. Further, 
we decompose Sjff into the orthogonal sum 

^N= ® (2-5) 

Telrrr 

where S)'^^ consists of functions form S^f^ which transform under according to some 
irreducible representation T of V. The decomposition (12.51) reduces 'Hfj. We denote the 
selfadjoint restrictions of Tif^ to by ■ The spectrum of T-L% is the union of the 
spectra of i^^'^, T G Irr F. 

Together with the whole system of A^ particles we will consider its decompositions 
into two clusters. Such decompositions play an important role in the characterization 
of the essential spectrum of the operators S)^^ . Let Z = [Zi, Z2) be a decomposition 
of the index set /:={!,..., A^} into two disjoint subsets: 

J = Zi u Z2, n Z2 = 0. 

Let 

Hz,! := 5^(/^n + K)+ Yl 

n<j 

nz,2 := ^ Z^n + J2 

n<j 

We introduce the operators corresponding to noninteracting clusters, 
cluster transferred far away from the sources of the external field: 

nz,j ■■= Azjnz,jAzj, in i^zj := ® A„L2(M^ C^), J = 1,2, (2.8) 



(2.6) 
(2.7) 

with the second 
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where 

Kz,j := TTA„ = ® An- 

For a given cluster decomposition Z = {Zi, Z2) we denote by P^^ and P^^ the 
projectors onto the irreducible representations Ej and Tj of the restrictions of 11 and F, 
respectively, to the cluster of particles indexed by Zj, j = 1,2. 

Given representations Tj and Ej, projector P'^iP^i = P^jp'^i reduces 'Hz,j- We 
denote the reduced operators in 

:= P^P^^^z,, 

by IH^z^f^ 1 ^^"^ define 

Kj {Z, Tj, Ej) := inf Spec H^'f' . (2.9) 

We write (Ti, Ei, T2, E2) -< (T, E) if the corresponding term cannot be omitted on 
the r. h. s. of 

(T2,_B2) 

without violation of the inclusion. For Z2 ^ let 

inf {x^{Z,T^,E,) + >i2{Z,T2,E2)}, Z^j^0, 

x(Z,T,P):= <;(^^'^^^^^'^^)^(^'^) (2.10) 
^MZ,T,E), ^1 = 0, 

and 

x{T, E) := mm{x{Z,T, E) : Z = {Z,, Z2), ^27^0}- (2.11) 

We are now ready to characterize the essential spectrum of in terms of cluster 
decompositions: 

Theorem 2.1. (Morozov [9], Theorem 6) For N ^ N let T be some irreducible 
representation ofT, and E some irreducible representation ofU, such that p'^p^ ^ 0. 
The essential spectrum ofTi^^ is [m:(T , E) , 00) . 

Thus the bottom of the essential spectrum is equal to the minimal energy which 
the system can have if some of the particles are transfered far away form other particles 
and sources of external field. We will omit the proof of the following simple proposition 
based on the positivity of the interaction potentials (12.31) . 

Proposition 2.2. It is enough to take the minimum in (12.111) over Z with Z2 = {n}, 
n = 1, . . . N . Moreover, for such Z , X2{Z, -, ■) in (I2.10p is equal to rrin, the mass of the 
particle in the second cluster. 

As shown in [6j, the Brown-Ravenhall operators can have eigenvalues below 
the essential spectrum. Note that in view of the decomposition (12.51) these eigenvalues 
can be embedded in the essential spectrum of 9)ff. 

Our main result is the following theorem. 
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Theorem 2.3. For N G N let T be some irreducible representation ofT, and E some 
irreducible representation ofli, such that P^P^ ^ 0. Let (p be an eigenfunction ofTt^^ 
corresponding to an eigenvalue A below the essential spectrum, i.e. 

n^^^(j) = X(f), X<x{T,E). 

Then there exists S > independent of A and such that for 

it holds 

f e'^l^lU(X)|'dX < oo. (2.12) 

Note that for A close to the bottom of the essential spectrum s behaves linearly 
in (^x{T,E) — a). However, for Schrodinger [10], Dirac [TU], Chandrasekhar [T3], and 
one-particle Brown-Ravenhall operators [T7] s can be chosen to be proportional to 
the square root of this distance. This suggests a conjecture that for the multiparticle 
operators we are considering the actual rate of decay might have this property as well. 
But the proof of such a conjecture is yet obscure even in view of [T7j, since that result is 
obtained by comparison to the decay rate of the eigenfunctions of Dirac operator, which 
are nonexistent in the multiparticle case. 

3. Some properties of the model 

In this section we single out some simple properties of the multiparticle Brown- 
Ravenhall operators introduced in the previous section. The reason for doing so is 
twofold. First, it will allow the reader to see which properties are required in each 
step of the subsequent proof of the exponential decay. Second, this will allow us to 
reformulate the main result without referring to the explicit form of the potentials (12.21) 
and ( 12. 3p . thus making future generalizations easier. 

We need a bit of notation. Let {Qj}jLi be a collection of uniformly C^-regular 

N 

domains in R with bounded boundaries. For = 1, . . . , A^, s G M, and = x I],- we 

j=i 

introduce the anisotropic Sobolev spaces 

H'^{n,c^ )■■={ ® L2{n,,c^)) ®H%n„.,c^) ® ( ® L2(^-,c^)). 

j=l j=n+l 

Property 3.1. For any R> there exists a finite Cr^ such that 




In other words, the interaction potentials are locally square integrable. 
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Property 3.2. The external field potentials decay at infinity in the Loo-norm: 

lim ess supl V"„(x) I = 0, n = l,...,N. (3.1) 

Property 3.3. For any e > there exists R > big enough such that for all 
n<j = l,...,N 

ll^ni^||L2(R3iVn{|x„-xj|>Rj) < e min ||^|| ^1/2.^^3^ ^4^., for all tjj E if^/^(M^^, C^^). 



Proof. This follows from the weaker property 

lim ess sup |f/„j(x)| =0, n, j = 1, . . . , N 

of the potentials (12.31) . □ 
Property 3.4. The interparticle interaction potentials are nonnegative: 

Unj > 0, for all n<j = l,...,N. (3.2) 

This follows from the assumption that all the particles of the system have electric charges 
of the same sign. 

Property 3.5. There exists C > such that for any n = 1, . . . , N 

\{Vn^,ij)\^Cy\\^^f,\m^.;2, for any y., ^ G ify2(M=^^, C^"), (3.3) 

and for any n < j = 1, . . . , N 

\{Un,^,ij)\^CMH./2\mm/2, for any e H'/\R'''' X'")- (3.4) 

Proof. Inequalities (13. 3p and (13.41) follow from Kato's inequality (see [20] and [21] . 
Theorem 2.9a.) □ 



Property 3.6. There exists C > such that for any n = 1,...,N and any ip G 

\\UnM\ ^ C min ||?/'||j:^i.^3^ ^^n.. (3.5) 



\\Ur,.M 

k=n,j "k'^ 

It is not surprising to have the minimum on the r. h. s. of (13.50 . since Unj only depends 
on the difference x„ — Xj. Note that (13.51) can be applied even if ip is only known to 
belong either to HI{M:^^ , C^^) or to Hj{R^^, C^"). 

Proof. Inequality (13.51) follows from Hardy's inequality (see e.g. [22], page 55) and 
the properties of symmetric-decreasing rearrangements (see e.g. [23], Lemma 7.17 and 
relation (3.3.4)). □ 
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Property 3.7. There exist Ci > and C2 G M such that for any cluster decomposition 
Z 



for any iP e ® A„f^^/2(M^ C^), j = l,2. 



Proof. This is where we need the subcriticahty condition (12.41) . According to the result 
of [1], inequahty (12.41) imphes (13. 7p if = 1. For A^ > 1 it is enough to use Proprety 13.41 
to estimate 7iz,j from below by a direct sum of one-particle operators. □ 



Remark 3.8. By Properties \3.5\ and \3.7\ the quadratic forms of operators (12.81) (and, 



in particular, Hn) CLre bounded below and closed on (g) A„iy^/^(R'^, C^). Thus these 

operators are well-defined in the form sense, and so are their restrictions to invariant 
subspaces. 

Remark 3.9. Suppose that the potentials ( \2.2\\ and (12.31) are replaced by operators of 
multiplication by some measurable hermitian matrix-valued functions such that K„ are 
the operators of multiplication of spin or coordinates ofn^^ particle by Ax A matrix-valued 
functions V^(x„), n = 1,...,N, and Unj are the operators of multiplication of spin or 
coordinates of n^^ and j*^ particles by 16 x 16 matrix-valued functions Unj{^n — ^j), 
n < j = 1, . . . , N . Then the statements of Theorems \2.1\ and \2.3\ remain valid provided 
Properties \3.1\ — \3. 1\ hold. Indeed, Properties \3.1\ — \3. 1\ imply Assumptions 1 — 5 of 



which form the hypothesis of Theorem 6 of And in the proof of Theorem \2.3\ we will 
not need the explicit expressions (12. 2p and (12.31) . but only the properties listed in this 
section. 

4. Proof of Theorem D 

Some constants in the proof can depend on the masses of the particles. Since we only 
deal with a finite number of particles with positive masses, such dependence will not be 
indicated explicitly. 

Lemma 4.1. Suppose that for some a > 

e2"l^"IU(X)PdX < cx), n = l,...,N. (4.1) 



Then fl212|) holds with s = N'^/'^a. 
Proof. 



N N 
n=l n=l 

Thus (14.11) implies (I2.12p after summation in n. □ 



Exponential decay of eigenfunctions of Brown-Ravenhall operators 9 

It remains to prove that (14.1 p holds with some suitable a > 0. Without loss of 
generality we will consider the case n = 1. 
Let p G C^([0, oo), [0, oo)) be given by 



r ^3 



p{z) := { 



1 

2;--, zG[l,oo). 



(4.2) 



For e > let 



Note that for any e > 



/(X) := /(xO := (4.3) 
l + ep(|xi|) 



l|V/|U^ < 1. (4.4) 
Since G L2(M3^, C^'^), for n = 1 (jH]) is equivalent to 

N'''^<^llL2(IR3iv,c4'^) ^ (4.5) 

with C independent of e. Note that for any e > the function e"-^ is twice differentiable 
with bounded derivatives. Hence multiplication by e"-^ is a bounded operator in the 
Sobolev spaces if"(M^C^) with s G [0,2]. 

The following two lemmata will be important in the subsequent proof. 

Lemma 4.2. For any G [0, 1) there exists C(ao) > such that for any a G [0, Oq] 
and%l)& L2(M^C^) 

||[Ai,e'^^]^L,(^3,c4)^C'(ao)a||e'^^^||, (4.6) 

and 

||e-'^^[Ai,e'^^]^||^,(^3,c:4) ^ C'(ao)a||^||. (4.7) 

Lemma 14^2] is proved in Section [6l Some analogous estimates with L2-norms instead 
of if^-norms can be found in [16]. 



Corollary 4.3. For any Oq G [0, 1) there exists C(ao) > such that for any a G [0, Oq] 

andi) G L2(M^C^) 

||e-"^Aie'^^7AKC(ao)||^||. (4.8) 

Proof. 

e-'^/A^e'^/ = Ai +e-'^^[Ai,e'^^], 
and (!47D implies (iSl). □ 

Lemma 4.4. Let he the hall of radius R > Q inM? centred at the origin. For any 
a G [0, 1/2) there exist C{R) > and C{a, R) > such that for any ip G H^/^{R^, C^) 

\\MHmmBn,c^) < C^(^)ll^llHi/^(iJ3..c4) + C^(a,i?)l|e-''^^^||L,(R3,c4). (4.9) 
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We prove Lemma 14.41 in Section [71 

In order to be able to apply Lemma [4.41 we will only consider a G [0, 1/2). We can 
thus fix ao € [1/2, 1) and no longer trace the dependence of the constants in Lemma [12] 
and Corollary 14.31 on this parameter. 

Let us fix a cluster decomposition 

Zo:= ({2,...,iV},{l}). (4.10) 

Then 

A^e"^0 = P^P^A^e'^/^ = ^ {P'^^ P^^ P^^)K^e''f (j). (4.11) 

(Ti,£;i;r2,l)-<(T,£;) 



The eigenfunction belongs to the form domain of Tijy' , which is 
P^P^ S Ar.H'/'iW'X') C H'/^iR'^'X''')- 

n=l 

Hence by ([TOD . (M . (IXTOD . and (l2ll]) 

(Aie'^^0, (Hzo,i+Hz„,2)Aie'^^0) 



^ (Aie'^^^, 5^ (xi(Zo,Ti,Pi) + X2(Zo,r2,l))(P^ip^i®P^^)Aie'^^^ 

Ei;T2,l)- 

2 



{TuEv,T2,l)-<{T,E) 
Zo 



(4.12) 



Let us introduce 



:= x(T,E)(e'^^<^, [e'^^ Ai]</.), (4.13) 

Q2:=(Aie"^0, (5^(/^n + K)+ 5^ t/„,+Di)[Ai,e"^]0), (4.14) 

n=2 l<n<j 

Q3:=(Aie"^0,[/^i,e'^^]0), (4.15) 

Q4 := - (Aie'^/Aie'^/0, (^i + J^t/i,)^). (4.16) 

i=2 
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Then by ( KWf (recall the definitions O, (EZD, dZBD, and flO) ) 
x(T,E)||e"^0f = (Aie'^^0,x(r,E)Aie'^^0) +gi 

AT N 

= (Aie'^^^, + K) + Yl ^-J- + ^O^"^*^) + Qi + Q2 

n=2 l<n<j 

Af AT 3 

= (Aie'^^0,e'^^(^(D„ + K)+ Yl Unj+Diy) + J2Qi (4-17) 

n=2 KrKj (=1 

4 4 

= (Aie'^^(^,e'^^H^'^0) + J2Q^ = M\^i^^^<Pf + 

z=i z=i 

4 



Z=l 



Thus 



(x(T,E)-A)||e»/0||H5^Q,, (4.18) 

1=1 

and it remains to estimate Qi, . . . , Q4. This will be done in the next four lemmata. 
Lemma 4.5. There exists a positive constant Ci such that 

IQil ^Cia||e'^/0f. (4.19) 

Proof. By (14. 13^ and Lemma 14.21 we have 

IQil ^ |x(T,E)|||e'^^0||||[e'^^Ai]0|| ^ C a\x{T , E)\\\e''^\\\ 

□ 

Lemma 4.6. There exists a positive constant C2 such that 

m^C^aWe'^^cPf. (4.20) 

Proof. Since Ai commutes with J2n=2i^n + ^n) + Z]i^n<j ^"i' ^ = 
Ai[Ai,e''-'']Ai = 0, we have 

. N Af X 

(Aie'^^0, ( 5^(D„ + K) + ^ f/„, j [Ai, e-^^l^) = 0. (4.21) 

^ n=2 l<n<j ' 

According to Lemma 14.21 

|(Aie"^0,Di[Ai,e»^]0)| ^ ||Aie"^0|| || [Ai, e»^]0|| ^ Ca||e"^0f . 
By dUD and (H:2T|) this implies fl4:20D . □ 
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Lemma 4.7. There exists a positive constant C3 such that 

IQ3I ^C3a||e"^0f. (4.22) 

Proof We have [Di, e^f] = [-ia ■ V, e"-^] = -ia ■ (Ve"^) = -ia ■ a(V/)e"^. Now (021) 
follows from (HlHIl and fj44ll . □ 

Lemma 4.8. There exist C4 > and Co(a) > such that 

Q4 ^ C4a||e'^^0f + Co(a)||0||^,/,. (4.23) 

We give a proof of Lemma 14.81 in Section [HI 

Substituting the estimates (I419D . (ICTD . (|4:22D . and (14:231) into KTE^ . we conclude 

that 

4 

(x(T,E) - A - a5^Q)||e"^0f ^ Co(a) (4.24) 

Now if 



a < min I i ( 5^ q) (x(r, i?) - A) | , 



then the expression in brackets on the 1. h. s. of (14.241) is positive, and (14.241) implies 
(14. 5 p with a finite C independent of e. Theorem 12.31 is proved. 

5. Boundedness of integral operators 

In this section we collect some auxiliary material for the subsequent proofs of 
Lemmata 14.21 14.41 and 14.81 In order to be able to obtain the information on the 
boundedness of (singular) integral operators we will need the following two theorems: 

Theorem 5.1. (Stein [24j, Chapter 2, Section 3.2) Let K : C be a measurable 

function such that for some B > 

|ir(x)| ^ fi|x|-", |V/r(x)| ^ B\x\-'^-\ for almost every x G M", 

and 

I i\:(x)d"x = 0, for all < Ri < R2 < 00. 

Jri<\x\<R2 

For g e Lp{W), 1 <p< 00, let 

A(^)(x) := / i^(x - y)^?(y)d"y, e > 0. 

J|x-y|>e 

Then 

\\A{g)\l^BM\p (5.1) 

with Bp independent of g and e. 
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Remark 5.2. Inequality fl5.ll) shows that the operator A := lim exists as a bounded 

e—f+O 

operator in Lp(M'^) and its norm satisfies \\A\\p ^ Bp. 
The second theorem is known as Schur's test: 

Theorem 5.3. Let and {Q2, ^^2) be two spaces with measures. Let A{-,-) be a 

measurable (matrix) function on Qi x ^2 satisfying 

Ml := sup / |A(x, y)|d;Ui(x) < cx), M2 := sup / | A(x, y) |d/i2(y) < 00. 
Then the integral operator 

(A^)(x):= / A(x,y)^(y)d/i2(y) 
Jn2 



is bounded from L2{^2) to L2(f2i) and \\A\\ ^ -y/ Mi M2 . 

We will only use Theorem 15.31 in the case Vti = 0.2 = with Lebesgue measure. 

Note that in the case of convolution (i.e. for y4(x, y) = A(x — y), ^21 = ^2 = M*^) 
Theorem 15.31 reduces to Young's inequality for convolution with Li-function (see e. g. 
[25]). 

For a 4 X 4 measurable matrix function A on X W we define the corresponding 
integral operator by 

(A^)(x) := hm / A(x,y)^?(y)dy, g&Cl{R\C'). (5.2) 

We will only work with such A for which (15.21) is well defined and extends to a bounded 
operator in L2(M^, C^) either by Theorem 15. II (in which case A(x, y) has to depend only 
on (x — y)), or by Theorem 15.31 

In particular, according to the definition given above and Appendix B of |6j, the 
integral kernel of (A^ — 1/2) is 



;C(x, y) = /C(x - y) := —^-—^Kr (m|x - y |) 



+ Plf r^V''^ 1'" + ^T^^V^ ^oHx - y|) 
47r^ V |x — y I |x — y|^ ^ 



Ki(m|x — y|) zq;-(x — y) 



(5.3) 



The boundedness follows from Theorem 15.11 and (1A.2P . 

Note that the function (15.31) rapidly decays together with its derivatives if |x — y| 
becomes big. Namely, if for r > we define 



G(r):= sup |/C(x,y)|+ sup |Vx/C(x,y)|, (5.4) 

|x-y|>r |x-y|>r 

then by (IA.2p and the first asymptotic in flA.ip . for any i? > there exists C{R) > 
such that 

G{r) ^ C(i?)r-3/2e-'-, for all r ^ R. (5.5) 
We will also use the following elementary lemma (Lemma 10 of [9J): 
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Lemma 5.4. For any d,k & N there exists C > such that for any bounded 
differentiable function x on with bounded gradient and u E H^^'^(M.^,C'') 

IIXw||//i/2{Md,Cfe) ^ C{\\x\\L^(Rd) + ||Vx||L^(Kd))||M||^i/2(Kd^C'=)- 

6. Proof of Lemma 14.21 

To prove (14. 6 p it is enough to show that [Ai, e"-^]e~"-'^ is a bounded operator from 
L2(M^C^) to H\R^X^) satisfying 

l|[^i'e"^]e""lL,(M3,c4HHi(R3,c4) ^ Ciao)a, a G [0, 1). (6.1) 
The integral kernel of [Ai, e"-^]e""^ = [(Ai - 1/2), e^^je""^ is given by (see (I5l3l) ) 

([Ai,e"^]e-'^0(^'y) = /C(x,y)(l -e"(^W-^(^))), (6.2) 
and its gradient in x is 

(V[Ai,e'^^]e-'^0(^'y) = (V./C)(x,y)(l -e'^^^W-^^^))) 
+ a/C(x, y ) (1 - e'^^^ W-^^'^)) ) ( V/) (x) - a/C (x, y ) ( V/) (x) . 

We rewrite 

1 _ ea(/(x)-/(y)) ^ ■ (x - y) + i?i(x, y) + i?2(x, y), (6.4) 

where 

i?i(x,y) := l + a(/(x) -/(y)) _ e'^(/W-/(y)), 

and 

i?2(x, y) := a((V/)(y) ■ (x - y) + /(y) - /(x)). 

Since 

|e^-l-2| ^ {e-2)z^ for \z\ ^ 1, 

by (14. 4p we have 

|i?i(x,y)| ^ (e-2)a2(/(x)-/(y))' ^ (e - 2)a2|x - yp, for |x-y| ^a-^ (6.5) 

On the other hand, since a < < 1, for |x — y| > a~5 the functions 

|/C(x,y)i?i(x,y)| and | Vx/C(x, y)i?i(x, y) | 

are integrable in x or y with the integrals bounded by C{ao)a, as follows from (15. 4p . 
(15. 5p . and (14. 4p . Since / G C^(M^), by the Taylor formula we have 

/(x) - /(y) = (V/)(y) ■ (x - y) + {{Vf) (^x + (1 - Oy) (x - y), (x - y))^^, 
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where Vf is the Hessian matrix (i. e. the matrix of the second partial derivatives of /) 
and ^ G [0, 1]. Hence 

|i?2(x,y)| =a ((P/)(ex+(l-0y)(x-y),(x-y))M3 ^ a\\VfUJ^ - y\\ (6.6) 



where is bounded uniformly in e by (14.31) and (14. 2p . Substituting (16.41) into 

(16. 2p and (16. 3p . and using the estimates (16. 5p — (16. 6p we obtain (16. ip by Theorems 15.11 
and 15. 3[ This completes the proof of (14. 6 p . 

The proof of (14. 7p is completely analogous since the integral kernel of 

e-"^[Ai,e"^] = e-''^[(Ai - 1/2), e"^] 

is 

/C(x,y)(e"(^(y)-^W) - 1) 

(compare with (16. 2p ). 



7. Proof of Lemma 14.41 

Let ri e C^{R^, [0,1]) with 



r/(x) 



0, X e B2R, 

1, xeR=^\fi3/?. 



Since Ai is a bounded operator in if^/^(M^, C'^), by Lemma [5.41 we have 

l|Ai^||HV2(ij^,c4) ^ ||Ai(l -r7)^||^i/2^^^^^4^ + \\AivHm/^(B^,,c^) 

By (15. 4p we can estimate the second term on the r. h. s. of (]7.ip as 



|A,-"^ 



^(x,y)??(y)^/'(y)dy 



Vx^(x,y)r/(y)^(y)dy 



(7.1) 



dx 



Br \y\>2R 
4 . 

^ -ttR sup 

3 xeB/i 



|y|>2i? 



^(x,y)r?(y)^(y)dy 



|y|>2R 



+ 



Vx^(x,y)r/(y)7/'(y)dy 



\y\>2R 



3 



sup \K{x,y)\ + sup |Vxi^(x,y)| ) |V'(y)|dy 

|y|>2R ^xeBfl -cR„ / 



f G{\y\-R)\i;{y)\dy 

\J\y\>2R 



^^R'( [ G'-'-{\y\-R)dy)( [ G'^''^{\y\ - R)\^{y)\'dy] . 

<3 \J|y|>2i? / V^|y|>2R / 
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Since a < 1/2 and /(x) ^ |x|, we conclude from (15.51) that there exists C{a,R) such 
that 

\\MvHHHBn,C*) ^ C(a,i?)||e-2'^^^||i,(K3,c4), 

and dSD follows by fTO) . 

8. Proof of Lemma 14.81 

For j = 2, . . . , N we have 



(Aie'^^Aie'^^</),f/i,0) = (f/i,-e»^(^, e'^^c^) 

+ (?7i,e-'^^[Ai,e'^^]Aie'^^0,e'^^0) + (f/i,[Ai, e'^^]0, e-^^ 



•1) 



The first term on the r. h. s. of (18.11) is nonnegative by (13.21) . Applying (13.51) . Lemma 1121 
and Schwarz inequality we can estimate the last two terms by Ca||e"'^0||^. Hence by 

drmD 

Q4 ^ Calle-^^^f + I (Aie'^^Aie'^^0, | (8.2) 

and it remains to estimate the last term on the r. h. s. of (18.21) . 

Let xi ^ (M^, [0, 1]) be a function supported in \ Bi such that it is equal to 
1 on \ B2. For > 1 let 

Xr(X) := Xr(xi) := Xi{^i/R)- 

We have 

I (Aie'^^Aie'^^^, V,<j)) \ ^ \ (e-'^^Aie''^Aie'^^0, XRVie^^<p) \ 
+ |((l-Xii)Aie"^Aie''^0,ri0)|. 

By Corollary 14. 3[ 

||e-"^Aie"^Aie"^(/)|| ^ C\\e^^(j)\\. (8.4) 
Since xr is supported outside Bj^, by (13.11) we have 

^ eiRW^H, e{R) 0. (8.5) 

According to (13. 3p . 

|((1 - XR)Me''^Aie^f<P,V,<P)\ ^ C\\{1 - XK)Aie'^^Aie"^0||^i/2 (8.6) 

Since (1 — xr) is a smooth function supported in {|xi| ^ 2i?}, by Lemmata 15.41 and 14.41 
we have 

11(1 - x/?)Aie»^Aie'^^0||^v. ^ C(i?)||Aie'^^Aie»^0||^V2(^^^,j,3iv-3,c4-) 
^ C{RWfA,e'^f^^^^.^^^^^^,^_,^^,.^ + C{a,R)\\e-'^fA,^^^ ^^'^^ 
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By Corollary 14.31 the second term on the r. h. s. of (18.71) can be estimated by C(a, -R)||0||. 
Applying Lemma 14.41 to the first term we obtain 

C(i?)||e'^/Aie"/0||^./.(^^^^^3.-3,c.-) 

^ C^(a,i?)l|Aie"^0||^i/2(^^^^jj3iv-3,c4-) 

f - f \°-°) 

^ C(a,i?)||e"-^0||^i/2(5^^^xK3iv-3,c4^) +C'(«,^)l|e <P\\ l^,(m.^n ^^^n ^ 

^ C(a,i?)||0||^l/2(]R3iV_C4^)- 

Thus by (EH) — (IHISD 

|((l-Xij)Aie'^^Aie'^^0,V^i0)| <C(a,i?)||0||^,/2. (8.9) 

Estimating the r. h. s. of (18.31) according to (18.41) . (18.51) . and (18. 9p and substituting the 
result into (18.21) we obtain 

Q4 ^ CaWe'^fcPW + Ce{R)\\e^fcpf + C(a, R)ml,/,. 

Choosing R so that s{R) ^ a we arrive at ( 14.23^ . Lemma [4.81 is proved. 
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Appendix A. Some properties of modified Bessel functions 

The modified Bessel (McDonald) functions are related to the Hankel functions by the 
formula 

These functions are positive and decreasing for z G (0, oo). Their asymptotics are (see 
[26] 8.446, 8.447.3, 8.451.6) 



(A.i; 



ir.(.)^y^e-(l + 0(l)), .-.+00; 

i^o(^) = - log ^ (1 + 0(1)), iri(^) = -(1 + 0(1)), 2^+0. 

z 

The derivatives of these functions are (see [26j 8.486.12, 8.486.18) 

K'^{z) = -Kr{z), K[{z) = -Ko{z) - -Ki{z), z G (0, 00). (A.2) 

2/ 
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